extern doublereal dnrm2_(integer *n, doublereal *x, integer *incx);
